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Abstract 


large space structural systems, due to their inherent flexibility 
and low mass to area ratio, are represented by large dimensional mathe- 
matical models. For implementation of the control laws for such systems 
a finite amount of time is required to evaluate the control signals; and 
this time delay may cause instability in the closed loop control system 
that was previously designed without taking the input delay into con- 
sideration. The stability analysis of a simple harmonic oscillator 
representing the equation of a single mode as a function of delay time 
is analyzed analytically and verified numerically. The effect of in- 
herent damping on the delay is also analyzed. The control problem with 
delayed input is also formulated in the discrete time domain. 

I. Introduction 

Large flexible space structures have been proposed for possible use 
in communications , electronic orbital based mail systems, and solar 
energy collection. ' 2 The size and the lew mass to area ratio of such 
systems warrant the consideration of the flexibility as the main contri- 
bution to the dynamics and control problem as compared to the inherently 
rigid nature of earlier spacecraft systems. For such large flexible sys- 
tems, both orientation and surface shape control may often be required. 

The equations of motion describing the shape of any large space 
structure are either represented by a few partial differential equations 
or a large number of ordinary differential equations. As the partial 
differential equations are difficult to solve for control system design 
purposes, the structural dynamics are corrmonly described using Finite 
Element Methods (FEM) . Two typical large space structures namely the 
Hoop/Column antenna^ and the Space Station initial operational configura- 
tion (IOC) 4 are both described using 672 degrees of freedom. Thus the 
dynamics of a large space structure can be written as^ : 
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where 


M = NXN mass/inertia synmetric matrix 
K = NXN stiffness synmetric matrix 

Z = NXl generailize coordinates representing the degrees 
of freedom 


u c~ 


influence of the external forces in each degree of 
freedom = B'U. 


With the modal transformation 

Z = <t> q 

and the properties of the modal transformation such as 


<J) T M4> = I 

iTvj, — 


2 2 


- h 


<j> K<{) = diag [u^ n _ 

and neglecting the higher modes, equation (1) can be written in standard 
state space form as 

X = AX + BU 


( 2 ) 


where 


X = 2nxl state vector representing modal coordinates 
and their velocities [q,q] T 


U = mxl 


A = 


0 


control vector 


II 

I nxn 


-Wi 


i 0 


2> 

11 ! 

n| 


system matrix 


B = 


0 


nxm 


4> T B * 


nxm 


control influence matrix 


II. Control with Delayed Input 

The proposed control systems for large space structures are based 
on state variable feedback of the form: 


U = -FX 


( 3 ) 
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and the control gain matrix, F, is designed using techniques such as the 
linear quadratic regulator (LQR) theory^, pole placement', and/or linear 
quadratic Gaussian/loop transfer recovery (LQG/LTR).® 


For the case when the^complete state is not available for feedback, 
an estimate of the state, X, is obtained using an appropriate estimator 
from the measurements of the form 


Y = CX 


(4) 


where 

Y = £xl measurement vector 
C = £xn sensor influence matrix 

In general, it is assumed that the estimated' state, X, is instanta- 
neously available. As the state estimator is implemented using a 
digital computer and the number of the status (2n) is of the order of 
hundreds for a large space structure, the computational time becomes 
appreciable. Thus, in the present paper, the stability of the closed 
loop control system, with the control as given in equation (2), is 
analyzed as a function of the delay time (h) using the modified control 
law of the form: 

U(t) = -FX(t-h) (5) 

The characteristic equation of the closed loop system 

X = AX(t) -BFX(t-h) (6) 

is given by 

— sh 

G(s,h) = det (sI-A+BFe ) = 0 (7) 

which, in turn, can be written as 

G(s,h) = P i (s)e“ shl = 0. (8) 

i=0 


The roots of the characteristic equation, (8), as a function of the delay, 
h, are obtained from the corresponding auxiliary equation^ 


where 


G' (s,h) 


£ P. (s) (1-Ts) 21 (1+Ts) 4n ' 2i = 0 
i=0 


-sh 

e 


[ 


1-sT 2 
1+sT J 


(9) 

( 10 ) 
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The value of T for which the roots of the equation (9) cross the 
imaginary axis in the s-plane is obtained and the corresponding h is 
evaluated using the relation, (10) . 

III. Example of a Harmonic Oscillator 


The equation of motion representing the ith structural mode is 
the familiar harmonic oscillator and is given by 


x i + “i x i - f i 


(id 


Considering the delayed velocity feedback of the form 


= “2£jU)^x^ (t-h) (12) 

with 

“i = 6 ' ?i = °- 5 ' 

the characteristic equation is given by 

G(s,h) = s 2 +36+6se -sh = 0 (13) 

= ilo p i( s ) e-shi = 0 

where 

P Q (s) = s 2 + 36 
P.,(s) = 6s 

The corresponding auxilary equation is given by 


E P. (S) (1-Ts) 2l (1+TS) 2-21 =0 (14) 

i=0 

(s 2 +36) (1 +Ts) 2 +6s (1-Ts) 2 = 0 

T 2 s 4 + (2T + 6T 2 ) s 3 + (l+ieT 2 -^) s 2 

+ (72T+6) s+36 =0 ( 15 ) 

Using the Routh-Hurwitz criterion, it can be found that the roots of 
equation (15) cross the imaginary axis at ui = 9.7 for T = 0.0426. 

The corresponding delay (h) can be calculated from the relation (10) 
with s = jw and is 0.16. This result can also be verified directly 
for this siirple system with the substitution s=jco into equation (13) 10 , 
resulting in the value of u> and delay h for which the roots of the 
characteristic equation cross the imaginary axis. 
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Thus, equation (13) can be written as (keeping and u> ) : 
(ok 2 -w 2 ) + j(2c i o) i U))e = 0 


(16) 


or 


2 2 

(u) . -a) +2c.u),o) sinwh) 4 * j2^.oj.aj cosu)h = 0 . 

l l i J i l 

For equation (17) to be satisfied 


(17) 


and 


or 


cos ooh = 0 or ojh = 


2 2 

a) . -w +2£.m.a) = 0 

i li 


(18) 

(19) 


Ol) - C.U). + 0) . / l+£ . 2 
11—1 1 


Taking the positive value for oj, the delay h, is given by 


tt I r 


h = 


j.[c.+ /TRT 2 " ] 

i i i J 


( 20 ) 


The value of h for = 0.5 and = 6 is 0.16 and thus the earlier 
result is verified. It is observed that an increase in damping reduces 
the tolerable delay (h) in the input. 

The equation of motion of a single mode with inherent (natural) 
damping and velocity feedback can be written as: 


X+2? 'ai .X+tu . 2 x = f = -lx, .O).x(t-h) 

111 11 

where is the inherent damping ratio. 

The corresponding characteristic equation is given by 

s +2£ !u) . s+u) . 2+2 £ .u.se = 0 . 
i 1 1 l i 

After substituting s = jw, equation (22) can be written as: 

2 2 ' 

(to . -co +2r, oj.uj sinoih) + i (2c .u ,u)+2£ .w .tocosuih) = 0 
i ii 1111 

For equation (23) to be satisfied for all u) and h, we have 

T 

2t .to .+2? .to . coswh = 0 
^i i ii 


or 


cos rnh = -c . / 
iC 


( 21 ) 


( 22 ) 


(23) 


(24) 

(25) 
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Thus, for coswh = < 1, the inherent damping must be less than damping due 
to control for instability. For the system will always be stable. 

With the value of uh from equation (25) the frequency cu can be 
calculated as: 


uj = oj [/ £ 2_£*2 + / I+ 5 2 — 5' 2 3 (26) 

i i i i 

and selecting the positive value of u), h is given by: 
cos 

h = ; 


u [✓ &'2 + ✓ 1+Uc ' 2 (27) 

1 i i i i 

For C. = 5 ' it can be seen that the delay, h, is half the undamped 
natural period of vibration. As the damping due to control increases, 
the tolerable delay (h) decreases and is in accordance with the observa- 
tion made in the case without the inherent damping. The effect of inher- 
ent damping in the system is to increase the amount of delay that the 
system can tolerate without become unstable as compared to the case 
without damping 


IV. Discrete Time Domain 


As the controller is implemented on a digital computer, it may be 
more natural to consider the delayed input problem in the discrete time 
domain . 

The equations of motion as given by equation (2) can be written in 
the discrete time domain as 


X(i+1) = A d X(i)+B d U(i) 


where 


A = e AA , B - / e A( ' t ' A) B dt 

d d 

o 

A = discretization time. 


(28) 


The delayed input . problem can be considered in discrete time in one of 
the two following ways: 

i) Designing the controller of the form U = -FX(i) without taking 
into consideration the delay and then examining the effect 
of delay on the stability of the closed-loop control system. 

The control gain matrix F is designed such that the matrix (A^-B^F) 
has the eigenvalues within the unit circle. Then the delay is introduced 
into the control law as: 
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U(i) = - FX(i-Jl) 


(29) 


and 


X(i+1) = A d X(i)-B d FX(i-£) . 


(30) 


The stability of equation (30) can be studied using the augmented system 
given by 


x(i+l) 


X(i) 

= 

X(i-S,+1) 



Z(i+1) 


A d 0 0 0 -B d F 

1 0 0 0 0 


0 0 0 


Oi 

A , 


X(i) 

X(i-l) 

X(i-Jl) 


% 

Z(i) 


(31) 


or 


'U 


Z (i+1) = A d z(i) (25) 

(ii) Designing the control by taking into account the delay in 
the input. 6 ’ 11 

Equation (28) can be modified as : 

X(i+1) = A d X(i) + B d U(i-i) 

The control law of the form U(i) = -FZ(i) can be designed from the 
augmented system: 


(32) 


X(i+1) 

U(i) 

U(i-l) 

U(i-£+l) 


A, 0 0 

d 


0 0 
I 0 


0 0 0 


0 B 

c 

0 0 

0 0 


X(i) 

U(i-l) 


U 


U-O J 


% 

Z(i) 


+ 


0 

I 

0 


U(i) 


(33) 


or 


Z(i+1) = A , Z ( i ) + Bn(i) . 
a a 
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Thus the input U(i-£) is a function of the previous inputs, U(i-£-l), 
U(i-£-2) , . . . , and the previous states X(i-£) . Though this design can 
take delay into consideration, the sequence of the control signals: 
U(i-£), U(i-£-KL) , . . . must be generated at an interval of one step and, 
thus, the original delay problem is not completely solved. 


Conclusions 


The effect of delay in the input on the stability of the continuous 
time controller that is designed without taking this delay into consider- 
ation is presented. The closed-loop control system of a second order 
plant becomes unstable for a delay of 0.16 seconds, which is only 16 per- 
cent of its natural period of motion. It is also observed that even a 
small amount of inherent (natural) damping in the system can increase the 
amount of delay that can be tolerated without the system becoming unsta- 
ble. The delay problem is formulated in the discrete time domain and an 
analysis procedure is suggested. 
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